Attractive Forces Between Global Monopoles and Domain Walls 
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We study the interaction between stable monopoles and domain walls in a SO(3) x Z{2) scalar field 
theory. Numerical simulations reveal that there is an attractive force between the monopole and the 
wall, but that after the monopole and the wall collide, the monopole does not unwind. We present 
an analytic explanation for the origin of the attractive force, and conclude that this is a generic 
feature of monopole- wall interactions which does not depend on the detailed structure of the m odel. 
The existence of the attractive force supports the hypothesis of Dvali et al. (hep-ph/9710301), who 
proposed that monopoles can be "swept up" by domain walls, thereby alleviating or solving the 
monopole problem associated with phase transitions occurring after or in the absence of inflation. 

BROWN-HET-1207 



I. INTRODUCTION 

In a previous paper in collaboration with Andrew Sorn- 
borger, we investigated the interaction of global topolog- 
ically stable monopoles with embedded domain walls in 
a 5'0(3) X Z2 linear sigma model |]. We were motivated 
by Dvali et aZ.'s proposal ||] that monopoles could be 
swept up by domain walls which, if true, provides a po- 
tential solution for the monopole problem. In our previ- 
ous analysis the Z2 symmetry factor corresponding to the 
domain wall was a "semi-discrete" symmetry, i.e. part of 
the connected symmetry group. Hence, the walls corre- 
sponding to this Z2 were not stable; they were embedded 
walls |3[Q]. However, we confirmed that both monopoles 
and anti-monopoles are attracted to the domain wall and 
that after colliding with the wall they unwind before the 
wall itself decays. Furthermore, we confirmed Dvali et 
aL's conjecture that the monopole charge dissipates on 
the surface of the wall by studying the evolution of both 
the winding number , and the winding number density 

Pw 

The question immediately prompted by our previous 
results is whether the same situation will be found when 
global monopoles interact with topologically stable do- 
main walls. Recent work provides numerical evi- 
dence that the attraction-unwinding mechanism extends 
to models with both topologically stable local monopoles 
and domain walls. In this paper, we ask whether or not 
the attraction- unwinding mechanism carries over to mod- 
els with stable global monopoles and Z2 walls. 

In our previous work we argued that the attrac- 
tive force between the monopole and the wall would 
also be found in systems with general global and gauge 
monopoles and domain walls because the underlying 
physics is generic, and rests primarily on the minimiza- 
tion of the gradient energies of the field configurations. In 
this paper we present evidence supporting this conjecture 
by studying a G = 0(3) x Z2 linear sigma model spon- 



taneously broken to H = 0{2). The vacuum manifold 
is M = G /H — S"^ X Z2 and thus admits both topologi- 
cally stable domain walls and monopoles. Our numerical 
simulations confirm the existence of the attractive force 
between monopole and domain wall. However, in con- 
trast to the interaction between global monopoles and 
embedded global walls, we find that the monopole fails 
to unwind on the surface of the domain wall. Finally, we 
provide an analytic explanation of the results seen in the 
numerical simulations. 



II. THE MODEL 

Topologically stable monopoles are ensured ||^ when 
the second homotopy group of the vacuum manifold is 
nontrivial, i.e. Il2{A4) ^ 0. Similarly, stable domain 
walls require the existence of a nontrivial Hq^M). There- 
fore, to get both stable domain walls and monopoles one 
needs Z2 and S'^ factors to be present in the vacuum 
manifold. In this investigation we consider a field theory 
with a novel potential which admits just such a vacuum 
manifold. The symmetry breaking pattern is 



G = 0(3) X Z2^H = 0(2) . 
Therefore, the vacuum manifold is 



G _ Q(3) X Z2 
H SO{2) 



X Z2 



(1) 



(2) 



To realize a symmetry breaking pattern which yields 
both stable monopoles and domain walls, we use the 
following Lagrangian: 
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After September 1, 2000: Imperial College, London. 
After October 1, 2000: Columbia University, New York. 



where = (01,02,03,04) is a four component scalar 
field. 

This model has an exact Z2 symmetry corresponding 
to the discrete transformation 04 — > —04. Furthermore, 
the Lagrangian is invariant under an 0(3) rotation of the 
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fields (jiij 4'2, (t>3- Hence the Lagrangian has an 0(3) x Z2 
symmetry. The set of vacuum states is 



J.2 2 
94 = '72 



(4) 



and the vacuum manifold hence is S*-^ x Z2. 
Since 

li^iM) = Ii2{GlH) = Z xZ2^l 



the existence of stable monopoles is ensured. A monopole 
configuration is 



^1,02,03) = f{r,t){-, -, -)?7i 



(6) 



where r — x + y + z and where the profile function 
f{r,t) obeys f{r,t) = f{r) as r and f{r,t) = 
/(r) — > 1 as r — > cx). A reasonable ansatz for the profile 
function at the initial time is 



(7) 



where Tc is the core radius of the monopole which is de- 
termined by balancing gradient and potential energies, 
which yields 
Likewise, since 



(8) 



our Lagrangian also supports the existence of stable do- 
main walls. A domain wall is constructed by picking one 
spatial direction (e.g. the z direction), and setting up a 
kink configuration in the 04 field: 



t>A{z) = rj2 tanh ( — 
V w 



(9) 



and letting the other three fields be at a fixed point in 
the S"^ subsp ace of A4. Here, w is the thickness of the 
wah, ^2/A2/?y2. 

This model has several advantages compared to the 
model with three scalar fields studied in In particu- 
lar, it admits stable domain walls rather than unstable 
embedded walls. Moreover, it is much easier to set up 
a configuration corresponding to a monopole displaced 
a certain distance zq from a domain wall as one does 
not have to rely on a product ansatz to simulate the 
monopole / domain wall configuration. The 04 field in- 
dependently generates a pure domain wall while the other 
three fields will support the monopole solution without 
affecting its winding number. Specifically, we can start 
with the (j)4 field in the configuration (j^) and set up the 
other three fields in the configuration (||) centered not 
at {x,y,z) = (0,0,0), but at {x,y,z) — (0,0, zq)- The 
mixing of the winding number due to the product ansatz 
was a concern in our previous work and now we can 
address it explicitly. 



We also need to study the monopole winding number, 
Qw, in order to determine whether the monopole un- 
winds. We can compute both the winding number and 
the winding number density. The topological charge as- 
sociated with Il2{A4) (which in our case is Z) is the wind- 
ing number. This number quantifies how many times the 
field configuration (j){x) wraps around the vacuum man- 
ifold as X ranges over a sphere 5^ in coordinate space. 
Hence the winding is defined as the homotopy class of 



(5) the map 
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from coordinate space to the vacuum 



manifold, known as isospace: : S^p^^^ Sf^o- 

For our 0(3) theory, the winding number in a volume 
enclosed by the surface Sk is: 



(10) 



where the indices a,b,c run from 1 to 3. 

Since we consider the time evolution of the winding 
over the entire coordinate space, useful information is 
obtained from the winding number density. Using Stokes' 
theorem in Eq. |lO| by performing a total derivative on 
the surface flux in the integrand, the winding number 
becomes: 



N : 



1 

8^ 



d''xeabce''''d,rdJ'>dkr. 



(11) 



By visualizing the evolution of the integrand in Eq. ^ 
we obtain information about the topological charge den- 
sity over the whole space. In principle, we could track 
the charge by studying the surface flux alone, but to 
do this one has to know where the winding is in order 
to choose a correct surface of integration, or compute 
the surface integral over many different trial boundaries. 
Consequently, we have found the winding number density 
to be a more useful variable. 

The one drawback of using the winding density to track 
the evolution of the monopole winding is that the density 
is initially a delta function, since eabci^-''^di(j)°-dj(j)'^dk(t>'^ ~ 
d{r) and the delta function is ill defined after we dis- 
cretize the equations of motion. Consequently, while the 
observed peak in the winding density will give an accu- 
rate qualitative picture of the monopole's location, the 
volume integral Eq. ^ is not conserved during the nu- 
merical evolution, even when we are looking at an iso- 
lated monopole. 



III. RESULTS AND DISCUSSION 

Partly as a test of the code, we first investigated the 
monopole and domain wall in isolation. In both cases we 
found that the solitons are stable. Note that if the initial 
thickness of the solitons is not chosen ideally, there will 
be excess energy which can and will radiate away. Due 
to the reflective boundary conditions implemented, this 
energy reflects back towards the core. Nonetheless, in 
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our studies of monopole - wall interactions this was not 
an obstruction, since the interaction occurred before the 
boundary effects became significant. 

Our primary goal was to investigate the attractive force 
between the monopole and wall. We did this by choos- 
ing initial conditions corresponding to a monopole that 
was initially at rest adjacent to a domain wall, and then 
numerically evolving the three dimensional equations of 
motion. The evolution of the energy and winding den- 
sities for a typical set of parameters is displayed in Figs 
1 and 2. The parameter choices for the plotted solution 
are Xi = rji = r]2 = 1 and A2 — 2, and the equations of 
motion were solved on a 192^ point grid. 

The attractive force between the monopole and wall 
is obvious from the simulation, and does not depend on 
any special parameter choices. In Fig. 1, we see that 
the monopole starts moving towards the domain wall, 
but does not pass through it. At later times (not shown 
in the plots) the distance between the monopole center 
and the wall does not decrease uniformly, but oscillates. 
This is apparently due to the natural oscillations of the 
monopole that are observed in isolated monopoles if the 
initial profile is not a time independent solution, and may 
possibly be accentuated by the edge effects due to the 
finite size of the simulation volume. However, there is 
no ambiguity about the existence of the attractive force 
that leads to the initial displacement of the monopole. 

We see from Fig. 2 that the monopole winding is ini- 
tially localized at the core of the monopole. The charge 
remains localized on the wall at the point of impact for 
late times and does not dissipate (unwind) on the wall. 
This differs from the evolution seen when a monopole 
interacts with an embedded wall |^], in which case the 
monopole can unwind on the wall. Consequently, we will 
now examine the differences between the two models so 
as to understand the lack of unwinding with a stable do- 
main wall, but in order to do this we must understand 
the mechanism which attracts the domain wall to the 
monopole. 

The mechanism underlying the attraction between gen- 
eral monopoles and domain walls is based on the assump- 
tion that the fields evolve to minimize the sum of gradi- 
ent and potential energies. As shown in [|], the gradient 
energy decreases as the monopole approaches the wall. 
Another way to understand the origin of the attraction 
is to focus on the change in potential energy and to asso- 
ciate it with a force between the monopole and the wall. 
If the force is attractive, it behaves like a tension which 
induces an acceleration along the line that connects the 
solitons. 

From the form of the potential in (|^) it is obvious that 
the potential energy associated with the monopole part 
of the configuration (0i — ^3) decreases with the dis- 
tance between the monopole and the domain wall, as the 
monopole's contribution to the potential energy being 
well localized. The potential energy is proportional to (/)| 
and thus decreases as the distance between the monopole 
and the wall decreases. However, this is a short range 




FIG. 1. Time evolution of the energy density (vertical axis) 
of a monopole colliding with a domain wall. The monopole is 
attracted towards the wall, and both the monopole and wall 
remain well defined after the initial interaction. The first 
frame shows the initial configuration, and subsequent frames 
are shown at times 1.42n (n = 1,2,3), measured in the same 
(arbitrary) units in which 771 = r;2 = 1. 
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FIG. 2. Time evolution of the winding number density of 
a monopole interacting with a domain walL The frames are 
synchronous with the energy density plots shown in the pre- 
vious figure. These plots show a detail from the center of the 
simulation volume, whereas the previous figure displays slices 
through the whole volume. The height of the peak is not sig- 
nificant, due to the inability of the discrete grid upon which 
the solutions are calculated to accurately depict an approxi- 
mate delta function. 



force since the monopole profile function /(r) approaches 
1 at large r. The decrease in the gradient energy as the 
distance between monopole and domain wall shrinks is a 
longer range effect (for global symmetries - for local sym- 
metries the gradients are screened by the gauge fields and 
hence also the force stemming from minimizing the gra- 
dient energy is short range). This effect comes from the 
fact that the volume of the region over which large gradi- 
ents are localized shrinks as the defects approach. Note 
that this argument works for general classes of defects, 
and adds support to the first crucial assumption in the 
DLV 1^ mechanism for solving the monopole problem. 

We can understand the lack of unwinding in the 
present model by studying the energetics associated with 
a monopole unwinding on the domain wall. Let us as- 
sume that the monopole winding can spread out to a 
distance I in the plane of the domain wall. When the 
monopole hits the wall its potential energy vanishes. 
What is left over is the gradient energy of both the 
monopole configuration (since the domain wall is topo- 
logical, there is no reason to expect the (/<4 configuration 
to change). As the monopole winding tries to spread out, 
the (jii — 4>3, hedgehog configuration will be distorted so 
that for most values of -\- < l"^ the 03 component 
dominates. Thus, in order to estimate the energetics as- 
sociated with a speading monopole charge, let us foucs 
on the contribution of the (/)3-component of the hedgehog 
field to the gradient energy. Making use of the symme- 
try in the plane of the wall, the energy E can thus be 
approximated by integrating the gradient energy density 
of the 03-component of the monopole integrated over all 
z out to a distance / in the plane of the wall. 

E > Ci^fdz ^ r -^dz (12) 
Jw oz z^n^ 

where w is the thickness of the wall along the z axis. 
Thus, 

E>-fri^w-^. (13) 

TT 

Hence we see that the energy grows as Z^. We can ob- 
tain an effective force by taking a derivative with respect 
to I. What we see is a restoring force F which acts to 
keep the monopole winding charge confined to a localized 
region of the domain wall. 

F = -kl with k = -rfw-^ . (14) 

TT 

We conclude that this "restoring force" is the basis of 
the non-unwinding of the monopole on the domain wall. 

IV. CONCLUSIONS 

We have studied the interaction between monopoles 
and domain walls in a scalar field theory with symme- 
try group G = 0(3) X Z2 symmetry group spontaneously 
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broken to = 0(2), a theory which admits both stable 
monopoles and domain waUs. We demonstrated the ex- 
istence of one of the two mechanisms which was crucial 
for the proposed solution of the monopole problem of , 
namely the attractive form between monopoles and do- 
main walls. However, in contrast to the case of global 
monopoles interacting with embedded walls ^ we found 
that the charge of the global monopole will not unwind on 
the surface of the domain wall. The gradient energy asso- 
ciated with delocalizing the monopole topological charge 
yields a tension which keeps the winding confined to a 
local region on the wall. Hence, we conclude that global 
monopoles can not unwind unless there are gauge fields 
to cancel the tension associated with the gradients of the 
scalar fields. Therefore, the lack of unwinding of global 
defects on domain walls is a universal behavior. 
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